Nonextensive local composition models in theories of solutions 
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Thermodynamic models present binary interaction parameters, based on the Boltzmann weight. 
Discrepancies from experimental data lead to empirically consider temperature dependence of the 
parameters, but these modifications keep unchanged the exponential nature of the equations. We 
replace the Boltzmann weight by the nonextensive Tsallis weight, and generalize three models for 
nonelectrolyte solutions that use the local composition hypothesis, namely Wilson's, NRTL, and 
UNIQUAC models. The proposed generalizations present a nonexponential dependence on the 
temperature, and relies on a theoretical basis of nonextensive statistical mechanics. The g-models 
present one extra binary parameter gy, that recover the original cases in the limit — > 1. Compar- 
ison with experimental data is illustrated with two examples of the activity coefficient of ethanol, 
infinitely diluted in toluene, and in decane. 

PACS numbers: 78.30cd, 82.60.Lf 



43 

Oh- 



43 

o 

O 

43 

Or 



o 
in 
o 

o 

CN 



X 



I. INTRODUCTION 

The basic hypothesis of local composition models, em- 
pirically introduced by Wilson [l| , and later used in other 
theories that followed, e.g. NRTL @, UNIQUAC 0, as- 
sumes that the composition in the vicinity of a central 
molecule differs from the bulk composition, and this local 
inhomogeneity strongly affects thermodynamic proper- 
ties of the solution. The description of binary and multi- 
component mixtures depends on interaction parameters, 
that are estimated from experimental data, and they usu- 
ally have the general form 



or even with a logarithmic term 



(i) 



with Aa,i j = ciij — a,jj , and is a molar potential energy 
of interaction between species i and j, with aji = a%j. 
For a comprehensive approach to theory of solutions, and 
fluid-phase equilibria in general, see Q. Aa^- is originally 
assumed to be constant. Extensions of the models relax 
this hypothesis, and consider Aa^ = Aay(T), according 
to various functions, for instance, a linear relation [2], 



Actij = Aoy ( o + Aa ijt i{T - T ), 
the inverse of the absolute temperature @, 

Adi,' i 

Aaij = Aa ijt0 H ^rS 

a combination of both linear and inverse relation 
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Different thermodynamic properties (heat of mixing, 
heat capacity, limiting activity coefficient, vapor-liquid 
equilibrium, liquid-liquid equilibrium, etc.), and/or dif- 
ferent chemical systems (presence of alcohols, hydro- 
gen bonding, etc.), may require different expressions for 
Aa^ {T) . The use of these equations also varies accord- 
ing to the considered model. Despite of the variety of 
functional forms, the exponential nature of the param- 
eter, that stems from the Boltzmann's weight, is kept 
unchanged by all models. Our proposal is to consider 
the nonextensive Tsallis weight, as a replacement for the 
Boltzmann weight. The departure from the exponential 
behavior is, thus, intrinsically originated from the dis- 
tribution, and not due to empirical modifications on the 
temperature dependence of the parameters. 

The paper is divided as follows: in Section|TI]we briefly 
present basic concepts of nonextensive statistical me- 
chanics, that will be used later. Section lllTl introduces the 
g-local composition model, with the nonextensive weight. 
Section ITVl applies the nonextensive local composition to 
Wilson's, NRTL and UNIQUAC models. Section |V] il- 
lustrates the effect of the nonextensive parameter q%j on 
the temperature, with instances of activity coefficient at 
infinite dilution. Finally, Section |VI] is dedicated to our 
conclusions and final remarks. 



II. BRIEF REMARKS ON NONEXTENSIVE 
STATISTICAL MECHANICS 

Despite the outstanding success of the Boltzmann- 
Gibbs statistical mechanics (BG), there are systems that 
are not properly described by the equations that emerge 
from this formalism. Along the last two decades there has 



2 



been continuously and increasingly developed the nonex- 
tensive statistical mechanics. Its starting point is the 
generalization of the concept of entropy (the Tsallis en- 
tropy) [|, 



1 



sr-W q 



(7) 



with pi the probability of the microscopic state i, W is 
the number of microscopic states, k is a positive constant, 
and q is the entropic index. If q — > 1, Eq. recovers 



the BG entropy, Sbg = Si 
S q is a generalization of Sbg- 



k'^2 i ^ In and thus 
Legendre transforms are 



preserved in nonextensive statistical mechanics [T(| . The 
g-entropy in the microcanonical ensemble (maximization 
of S q with equiprobabilities) Q is S q 
the g-logarithm is defined as 



fcln, W, 



where 



ln g x 



A-i _ i 



(8) 



The celebrated equation, engraved in Boltzmann's tomb- 
stone, S = klnW, is recovered at q —> 1. Maximization 
of S q , Eq. ([7]), with the constraint of constant generalized 
mean energy, leads to the canonical ensemble distribution 
for the energy (see [IS (3 f° r details), 



p(x) cc exp q (-/3 q Ei), 



(9) 



where j3 q is is the Lagrange parameter, that is related 
to the inverse temperature (/?i = l/(kT) in the BG 
formalism), Ei is the energy of the t-th state, and the 
g-exponential is precisely the inverse function of the q- 
logarithm, Eq. ©, 



exp x = [1 + (1 — q)x] 



(10) 



The symbol [A]+ stands for [A] + = A if A > 0, and 
= if A < 0. Equation © is the Tsallis weight, that 
is a generalization of the Boltzmann weight. The main 
difference between Tsallis and Boltzmann weights is that 
the former presents power law tails (long-lasting for q > 
1, and abruptly vanishing for q < 1), while the later has 
exponential tails. Equations © and (fTU|) present many 
similar properties of the usual logarithm and exponential, 
e.g. ln 9 1 = and exp ? = 1, Vg, and the derivative of 
the g-exponential is given by 



d(exp g x) 
dx 



(exp x) 



(11) 



In general, (exp^a;) 7^ exp q (ax), except for q = 1. The 
g-exponential and the g-logarithm functions lead to a 
nondistributive g-deformed algebra [Hj],[IiJ|. Generalized 
algebraic operations (g-addition x(S q y, g-difference xQ q y, 
q-product x <8> g y, g-ratio x 9 y) are defined as 



xQ q y\ 



y = x + y+(l- q)xy, 



l + (l-q)y 



(12) 



(13) 



x q y = [x 



= 1^-9 



-1]]-" (x,y>0), (14) 



y 
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111- (x,y>0). (15) 



With these g-operations, the g-exponential follows the 
properties: 



exp ? x exp g y 
exp ? x/ exp ? y 
exp 9 a; ® q exp g y 
exp 9 a; <Z> q exp g y 



exp g (x© g y), 
exp 9 (x Q q y), 
exp q (x + y), 
exp q (x - y). 



(16) 



The q-algebra has been applied in different contexts 
within nonextensive statistical mechanics. The g-product 
has been used in the generalizations of the central limit 
theorem and the Fourier transform [II] , [3 , [13 : an d 
there are evidences that it is connected to g-Gaussian 
distributions p(x) = A(g)v / ?exp g (— j3x 2 ) [lg|. Some 
properties of g-functions and g-algebra may be found at 
H3|;[13;[1I|'[23; anc l references therein. 

Nonextensive statistical mechanics is expected to be 
valid in a variety of situations: systems with long range 
interactions, long term memory, fractal structure, break 
of ergodicity, quasi-stationary states, or other features 
that characterize complex behavior. Let us briefly elabo- 
rate on the range of interactions, and on the nature of the 
quasi-stationary states, following the lines of [23j]. We can 
generally consider that interactions decay with distance 
r as l/r a . If a < 3, the interaction is long-ranged; a = 1, 
for Coulomb and gravitation interactions, are typical ex- 
amples. Such systems may exhibit negative specific heat, 
e.g. [24j . If a > 3, the interaction is short-ranged, e.g. 
van der Waals a = 6 interactions. Dipole interactions are 
at the threshold (a = 3). See [23[ for weak and strong 
violation of BG, and more details. 

Thermodynamic equilibrium is concerned about two 
limits: the time limit (time t — > 00, related to equilib- 
rium) and the macroscopic limit (number of particles 
N — > 00, related to thermodynamics). For simple sys- 
tems, the order in which these limits are taken is irrel- 
evant, but for complex systems, these limits may not 
commute. For certain classes of complex systems, if 
the time limit is taken first, and then the macroscopic 
limit (limAT^oo limt-^oo /(<, iV) of a dynamical function 
fit, N)), the system is characterized by Boltzmann equi- 
librium distributions, thus q = 1. But if the limits are 
taken in the reverse order (lim^oo limjv-foo f(t, N)), the 
system may achieve a quasi-stationary metaequilibrium 
state, according to its initial conditions, that is possibly 
described by nonextensive distributions. This hypothesis 
was conjectured in 1999 by Tsallis [23j, and it was compu- 
tationally verified for conservative long-range interacting 
systems, e.g. [H , [11 , [l?! . 

We address some examples that follow nonextensive 
behavior. The rate of re-association of CO with Myo- 
globin dissolved in glycerol-water solutions, after being 
photo-dissociated, was found to be described by func- 
tions that are connected to nonextensive statistical me- 
chanics [2||. This was supposed to be related to the 
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path CO molecules have to perform in a fractal-like struc- 
ture though the interior of folded proteins. g-Gaussian 
distributions were applied to polymeric networks, when 
finite chain effects are relevant 29]. Nonextensive dis- 
tributions of velocity of monomers during the relaxation 
process were found by molecular dynamics simulations 
of polymer chains and Lennard- Jones molecules [301 ] . q- 
Exponential functions were used to describe diatomic po- 
tential energy curves, particularly for and L12, and vi- 
brational spectra and spectroscopic constants were found 
to be in good agreement with experimental data (3~i| . Ar- 
rhcnius law has recently been generalized by the use of 
the ^-exponential, and agreement with experiment was 
found in plant respiration rates, bacterial gliding, and 
tunneling in the F + H2 reaction [32j. Theoretical and 
experimental aspects, and the historical development of 
nonextensive statistical thermodynamics, and many ex- 
amples, may be found in (22j. See also [33] for a theoret- 
ical treatment of generalized thermostatistics. 



III. NONEXTENSIVE LOCAL COMPOSITION 

The basic assumption of the local composition the- 
ory, as introduced by Wilson [l|, is that, due to dif- 
ferences in molecular sizes and in intermolecular inter- 
actions, the ratio of the number of molecules of species 
i and j, in the vicinity of a central molecule j, differs 
from that of the whole solution, according to Xij/xjj = 
(xi/xj)(exp(— exp(— ^)), where Xi is the bulk 
mole fraction of species i, and Xij is the local mole frac- 
tion of species i in the neighborhood of a molecule of 
species j. Our assumption is to replace the Boltzmann 
weight by the nonextensive Tsallis weight. Besides, the 
g-product of probabilities yields nonextensive distribu- 
tions (distributions remarkably close to g-Gaussians, to 
be more precise; see [l8[ , [34[ , [IH ) • The g-product, and 
the g-ratio, of Tsallis weights has recently been shown to 
be related to the locality of a generalized master equation 
35]. This inspires us to use the g-ratio, in the general- 
ization of the local composition hypothesis: for a central 
molecule j, 



with Aa,; 



cijj. Symmetry of the interactions 



Xij Xi { dij \ 

~ T = " cx P<?;i I — ~~EFr 



9iJ exp, 



'*■> V RTj 



(17) 



RTJ ^ 3 V RTJ 

Substitution of Eq. (|17p in the normalization condition 
12i x ij = 1 (c is the number of chemical species) results 

Xij = y^fr- (18) 

with the parameter A q> ij given by 1 

Aa,- 



An 



: exp 



RT 



(19) 



We adopt the symbol A q ^j, 
essary heavy notation. 



instead of A 



to avoid unnec- 



implies aji = ciij , and we assume, for simplicity, qji = qtj . 
The limiting case qtj — > 1 recovers the usual parameter 
Aij = Ai tij , Eq. ((J). Note that, according to Eq. (fTTj) , 



dA a 



dT 



A 



Aa* 



(20) 



Figures Q] and [2] show the behavior of the parameter 



A qi i2, Eq. (|19|) . as a function of Aai2, and the scaled in- 
verse temperature, respectively. The usual case q\2 = 1 
appears as straight lines in these semi-logarithmic plots, 
and q\2 > 1 (012 < 1) presents positive (negative) con- 
cavity. Equation © is also displayed in Fig. [1] for com- 
parison. The effect of the parameter Aai2,i of Eq. Q 
is simply to shift the q\2 = 1 curve, but it remains a 
straight line, once it relies on the Boltzmann weight. 
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FIG. 1: Parameter A q ^\2, Eq. (|19[) . as a function of the binary 
interaction parameter Aai2 (with T — 298.15 K). q±2 = 0.5 
(dot-dashed), 512 = 1 (solid), and qi2 = 1.5 (dashed). Dot- 
dot-dashed curve uses Eq. @ with Aai2,i = 5 x 10 5 J K/mol. 



IV. q- WILSON'S, g-NRTL, AND g-UNIQUAC 
MODELS 

A. q- Wilson's model 

The excess molar Gibbs free energy for the Wilson's 
model [l] is an empirical modification of the Flory- 
Huggins' equation (see [1] for details and additional ref- 
erences): 



9 

RT 



^7, 



(21) 



where t;u is the local volumetric fraction of the compo- 
nent i in the neighborhood of another molecule i (Flory- 
Huggins' model for athermal polymeric solutions uses the 
global segment fraction £j in the place of £u, in Eq. (|2"Tj) ). 
The local volumetric fraction is given by 



XaV-i 
j x ji v j 



(22) 



4 



10 



cr 10" 



10"' 



10' 



jr 10 r 





1 1 I 1 




(a) ; 




= 1.5 ■ 






! q = 0.\ 




■ Aa [2 = 5 kJ/mol 
, i.i. 


\ . i . 



2 3 4 

1 000/T (k 1 ) 





1 ' 


: Aa 12 = -5 kJ/mol / 


(b) ; 


- q=iy 




r / q ^ 




^7 1 .1.1, 


-q^05; 
i 



2 3 4 

1 000/T (k 1 ) 



10 L 



FIG. 2: Parameter v4 g ,i2, Eq. (|19|) . as a function of the inverse 
temperature (conveniently scaled), for two typical values of 
Aai 2 . The ordinary qi2 = 1 model appears as straight lines 
in these semi-log plots. q\2 > 1 {qi2 < 1) presents positive 
(negative) concavity. 



where the molar volume of liquid i, Uj, is taken as a 
measure of its molecular volume. 

Substitution of Eq.s ([18]) and j22]| in Eq. ([2T]l, leads to 
the nonextensive generalization of the Wilson's model, 
and the excess molar Gibbs free energy is 



g E c I c 



(23) 



with A gj y = (vj/vi)A q ,ij. The activity coefficient is 
formally given by the same expression of the original 
(q-ij = 1) model, but with the g-parameter A^-: 



In 7j = - In X! x 3 A q^3 + 1 - E 



(24) 



S/c X k^qjk 

The excess molar enthalpy (heat of mixing) h E = d Q 9 ^ 1 j^ 
is given by (see Eq. (|2TJ|) ) 

YT, -''/A,,. ,, 4^ V v. 



E 



El?) ^ A ^ A ^- ( 25 ) 



molar Gibbs free energy is given by 

c c 



(26) 



The NRTL model introduces a nonrandomness parame- 
ter a.ij — ctji , so the parameter A qi j is given by a varia- 



tion of Eq. (fll?|) . namely = exp^. . (— -^f). The 
expression for the activity coefficient of a component i is 
formally the same as the original (gy = 1) NRTL model, 
just replacing the usual parameter Ay by A q ^ . This pro- 
cedure is not valid to find the expressions of excess molar 
entropy and excess molar enthalpy, due to Eq. (|2H)) . 



C. q-UNIQUAC model 

The UNIQUAC model |[ is based on the local com- 
position hypothesis, but it replaces the local and global 
mole fractions in Eq.s (IT71) and (|18[) by the local and 
global surface fractions 0y and 0i (the global surface frac- 
tion is defined as 9i = Xiqij (Hj XjQj), where g, is the sur- 
face parameter 2 ). The interaction parameter ay is con- 
sidered as a measure of the internal energy (ay = ay), 
and the excess molar internal energy is then given by 

c c 

u E = Wi E ^fAtij-i. (27) 
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The excess molar Gibbs free energy is found by the ap- 
proximate relation 



9 E /T : 



1/T 



u E d(l/T), 



(28) 



that, for the multicomponent case, shall be numerically 
integrated. The residual contribution for the activity co- 
efficient of component i, for a multicomponent mixture, 
may be found by the numerical integration of 

ln 7 r s = C' T ufd{l/T), (29) 



9jAq,ij Auy 



with the partial molar excess internal energy given by 

c 

Qi E 



* ^ 121 9kA qM 



3 
c 

-Qi E ■ 

3 



ft A 

U 3^Q,V 



(121 °kA qtkj )' 



J2 c k k A q , kj Au kj . (30) 



B. g-NRTL model 

The nonextensive generalization of the NRTL model 
follows the same lines of [2j, with Eq. (|17p . The excess 



2 We use the notation <ji to avoid confusion with the nonextensive 
parameter qij. Sometimes there appear two surface parameters, 
q and q' , one of them is used in the combinatorial term and 
the other in the residual term of the UNIQUAC model pj. We 
consider q = q' , for simplicity. 
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Analytical solution for the binary case is as follows: where we have used the change of variables r = 

substitution of Eq. (EH) in Eq. Jgg), and Eq. QQ) in ex Pg (-cx), and f(r')dr' = ^ /(r')rfr' - f{r')dr', 

Eq. ([29]). with c = 2, lead to (see Eq. 3.194 5 and with 
Eq. 3.194 1 of 36]) 



cexp ? (— cx) 
o a + bexp q (—cx) 
1 



-dx 



2-9-, 



a{2-q) 

cexp (— car) 
o [a + &exp g (-cx)] 2 
1 



a Ji 1 
6 

— i 

a 



r' 1 


-Q 


1 + 


-7 C 




a 











dr' 



(31) 



« 2 (2-<7) 



dx = ;r 



(1 + V) 2 

"09 



X 9 (x) - 2^(1, 2 3- q;-x), (33) 



V>,(z) = 2F 1 (2,2-q;3-q;-x), (34) 



g < 2, and ^^a, /3; 7; x) is the hypergeometric func- 
tion, resulting the following analytical expressions for the 
,(32) residual contributions, g^ es and \njl es : 



RT 



Xiqi 



(2-? 13 )fii 

X 2 g2 #1 

'(2-? ia )fla 



2-q 12 10 2 \ I 2 

r q ,2i Xg 12 [ T r q ,2l J - X 9lm y ^ 



2-9l2 
T q,l2 Xqi 



-TqAI 



Xqi2 



(35) 



+ 



(2-9 la )»< 
_qi__0j_ 

(2-9 ia )«< 

(2-g 12 ) 



^912 I fl. T 9'J' i ) 



09 



Me 



(36) 



with g i2 < 2, r g ,y ee Aq,ij = exp q ..(-Auij/(RT)), 
and (i = 1, j = 2) or (i = 2, j = 1). The lim- 
iting case gi2 — > 1, with Xi(x) — x~ x ln(l + x), and 
tMx) = (1 + x)- 1 (see Eq^9.121 6 and Eq. 9.121 16 of 
[36j | . and Eq. 15.3.15 of [33]), recover the usual expres- 
sions. 

It is possible to find an expression for the residual con- 
tribution of the activity coefficient of species i, in a binary 
solution, different from, but equivalent to, Eq. (|36p . with 
the partial molar excess Gibbs free energy taken from 
Eq. ([35]). and hi7[ es = (gf) res /(RT), where it is neces- 
sary to use the derivative of Xq( x ) ( see Eq. 15.2.1 of 37]), 
d Xq (x)/dx = -(2 - g)/(3 - q) 2 F X (2, 3 - g; 4 - g; -x). 

The combinatorial contribution of the UNIQUAC 
model for the multicomponent case, that is the Guggen- 
heim expression for athermal mixtures, remains un- 
changed in the nonextensive generalization, consistent 
with the lower limit of Eq. ([28]) 0]: 

{Q com rP^ + lp^ (37) 



m7r m6 = ln <h + Z ln ^_ + k _^J2 Xj li (38) 
X% £ (pi Xi 

3 

where the volume fraction = XiTi/Q^ XjVj), is the 
volumetric parameter of molecule i, 1% = (z/2)(rj — qi) — 
(rj — 1), and z is the coordination number, usually as- 
sumed z = 10. The complete expression for the excess 
molar Gibbs free energy comprises the two contributions, 

9 9res ' 9 comb' 



V. ACTIVITY COEFFICIENT AT INFINITE 
DILUTION 

We have chosen the activity coefficient at infinite di- 
lution, that is a property that depends only on the tem- 
perature, to illustrate our proposal. We have also cho- 
sen the q- Wilson's model to illustrate the effect of the 
nonextensive parameter gy on the temperature depen- 
dence, because it is the simplest model, if compared 
to g-NRTL and g-UNIQUAC. The activity coefficient 
at infinite dilution for the g- Wilson's model is ln7f° = 



6 



In A g> i2 — A 9i 2i + 1. Fig. [3] shows curves for different val- 
ues of <7i2 (all curves use the same parameters Aa\2 and 
Ad2i). qi2 — 0.85 exhibits a maximum and a minimum 
in while q±2 = 1 can only display a maximum. 
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FIG. 3: 7i° according to q- Wilson's model, as a function of the 
inverse temperature, for typical values of the binary parame- 
ters: Aai2 = 10 4 J/mol and A<22i = —5600 J/mol (vi — V2, 
values of qi2 are indicated). Curve for qi2 = 0.85 displays 
a maximum and a minimum. The positions of the extreme 
points and the values of 7^ may be changed by adjusting 
the parameters to other values; the particular values of the 
parameters were chosen to make the figure visually good. 

We have fitted the q- Wilson's model to two examples 
of experimental activity coefficient at infinite dilution: 
ethanol infinitely diluted in toluene (Fig. |4^,) , and in de- 
cane (Fig. Hb). The examples are not properly described 
by the ordinary Wilson's model, and the generalized q- 
Wilson's model is able to describe the data. The data 
were taken from [H[ (Vol IX, Parts 3, 4) 3 . The two 
examples present positive deviations from Raoult's law 
(7i° > 1) and decreasing 7^ with the temperature. 



VI. FINAL REMARKS 

The main goal of this paper is to introduce nonex- 
tensivity in local composition models, used in theories 
of solutions. Current theories and models are based on 
Boltzmann's distribution and Boltzmann's weight. De- 
viations from Boltzmann's weight have been proposed, 
mainly on an empirical basis (see Eq.s (H])-©). The q- 
local composition hypothesis introduces one additional 
binary parameter that comes from the nonextensive the- 
ory, and generalizes the temperature dependence of the 
models. 



3 Four references are reported in [38ll for the system ethanol- 
toluene, and their fluctuations are significative. To avoid such 
fluctuations, we have considered one single set of experimental 
data, measured with the dilutor technique: Ref. 10 of [3S| . Vol 
IX, Part 3, pp. 1292-1293 
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FIG. 4: Activity coefficient at infinite dilution of ethanol in 
the solvents toluene (a), and decane (b). Experimental data 
from [3^. Dashed lines are the best fittings for the origi- 
nal (qi2 — 1) Wilson's model, and solid lines are the best 
fittings for the q- Wilson's model, (a) Dashed line: Aai2 = 
6880.0 J/mol, Aa 2 i = -582.0 J/mol; solid line: q 12 = 0.83, 
Aai2 = 10267.0 J/mol, Aa 2 i = -5179.7 J/mol. (b) Dashed 
line: Aai2 = 8773.6 J/mol, Aa2i = —26.6 J/mol; solid line: 
qi2 = 0.7, Aa i2 = 6500.0 J/mol, Aa 2 i = -2353.7 J/mol. 



As one increases the degree of freedom of a model by 
adding extra parameters, it is natural to expect a better 
optimization from a fitting procedure. Sometimes new 
parameters are simply introduced to take advantage of 
the additional degrees of freedom, and then turn the fit- 
tings easier. Nonextensive distributions do present a new 
parameter q^ , and of course the fittings are benefited 
from it. But the extra parameter qij lies on a theoretical 
background, and there are plenty of examples showing 
that the entropic index q has a physical interpretation, 
expressing the degree of nonextensivity of the system [22j . 

We have shown two fittings for the activity coefficient 
at infinite dilution for the q- Wilson's model. Expressions 
for other models (g-NRTL, q-UNIQUAC) were presented, 
but applications to large amounts of experimental data 
still remain to be explored. 

We hope this work invites experimentalists to apply 
nonextensive local composition models to large experi- 
mental databases, and test their validity in the descrip- 
tion, or predicting capability, of thermodynamic proper- 
ties, in different temperatures. 
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